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ABSTRACT 


Modern  astrogeodet ic  methods,  although  accurate,  are 
inefficient  and  too  complex  to  rapidly  determine  deflections  of 
the  vertical.  This  problem  is  the  impetus  for  finding  a  more 
useful  technique  that  can  yield  results  much  more  quickly.  The 
Global  Positioning  System  (GPS),  with  its  ability  tn  provide 
coordinate  differences  in  interferometric  modes,  can  be  used  to 
determine  these  deflections.  Using  highly  accurate  coordinate 
differences  in  conjunction  with  ortho me trie  height  differences, 
one  can  develop  a  surface  of  the  geoidal  undulations  as  a  function 
of  latitude  and  longitude.  Given  three  GPS  stations,  a  local 
surface  can  only  be  approximated  by  a  plane.  With  moie  points, 
however,  the  modeled  surface  will  more  accurately  resemble  the 
true  undulation  differences.  From  this  modeled  surface,  one  uses 
least  squares  fitting  of  polynomials  to  interpolate  the  "average" 
and  in  the  survey  area.  Finally,  the  deflections  and 

are  computed  and  a  study  of  propagation  of  both  absolute  and 
relative  errors  is  made. 
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ABSTRACT 


Modern  astrogeodet ic  methods,  although  accurate,  are 
inefficient  and  too  complex  to  rapidly  determine  deflections  of 
the  vertical.  This  problem  is  the  impetus  for  finding  a  more 
useful  technique  that  can  yield  results  much  more  quickly.  The 
Global  Positioning  System  (GPS),  with  its  ability  to  provide 
coordinate  differences  in  interferometric  modes,  can  be  used  to 
determine  these  deflections.  Using  highly  accurate  coordinate 
differences  in  conjunction  with  orthometric  height  differences, 
one  can  develop  a  surface  of  the  geoidal  undulations  as  a  function 
of  latitude  and  longitude.  Given  three  GPS  stations,  a  local 
surface  can  only  be  approximated  by  a  plane.  With  more  points, 
however,  the  modeled  surface  will  more  accurately  resemble  the 
true  undulation  differences.  From  this  modeled  surface,  one  uses 
least  squares  fitting  of  polynomials  to  interpolate  the  "average'1 
yy  and  y^  in  the  survey  area.  .  Finally,  the  deflections  f  and  v 
are  computed  and  a  study  of  propagation  of  both  absolute  and 
relative  errors  is  made. 
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INTRODUCTION 


Current  ast rogeodet i c  techniques  employ  ue 1 1 -known  geodetic 
concepts  uith  extremely  accurate,  complex  and  expensive 
astronomical  equipment.  Combining  geodetic  coordinates  from 
conventional  surveys  such  as  tr i angu 1  at i on ,  t r i 1 aterat i on  or 
traverse  uith  astronomically  observed  natural  coordinates,  one  can 
determine  deflections  of  the  vertical  for  many  points  on  the  solid 
surface  of  the  earth.  The  relationship  between  geodetic  and 
natural  coordinates  is: 

♦  =  ♦  -  {  (1) 

X  =  A  -  Tj  sec  ♦ 

(Heiskanen  and  Moritz,  1967)  where  (  is  the  deflection  in  the 
north-south  direction  and  v  is  the  deflection  in  the  east -west 
direction. 

Obtaining  accurate  geodetic  coordinates  uith  conventional 
methods  is  not  an  easy  task.  A  skilled  crew  of  four  to  five 
surveyors  must  spend  a  significant  amount  of  time  "in  the  field" 
connecting  previously-established  control  to  the  point  or  points 
of  interest.  Depending  on  the  distance  between  the  existing 
control  and  the  new  stations  and  on  the  roughness  and  vegetation 
of  the  terrain,  the  crew  may  need  many  days  or  even  weeks  to 


obtain  the  required  data.  Furthermore,  observing  stars  to 


determine  natural  coordinates  with  a s t r ogeode t 1 c  equipment  can  be 
accurately  completed  by  only  the  most  accomplished  of  astronomic 
surveyors.  Because  the  equipment  (usually  a  Wild  T-4  in  the 
United  States)  is  so  heavy  and  sensitive,  many  areas  do  not 
qualify  for  astronomic  observations.  Rough  mountainous  regions, 
swampland  and  glaciated  areas  cannot  adequately  provide  support  or 
shelter  for  the  telescopes.  Urban  areas,  polluted  skylines  or 
hazy  atmospheres  are  also  generally  out  of  the  astronomer’s  range 
because  of  the  inadequacy  of  the  light  coming  in  from  the  stars. 
All  in  all,  there  is  no  question  that  modern  astrogeodet ic 
techniques  are  burdened  by  a  plethora  of  difficulties. 

Observations  with  the  Global  Positioning  System  (GPS),  on  the 
other  hand,  are  tremendously  advantageous  in  that  they  defeat  many 
of  the  problems  that  plague  astronomic  techniques.  Firstly,  GPS 
receivers  are  continuously  being  ruggedized  and  can  already 
withstand  conditions  many  times  more  severe  than  what  an 
astr otheodo 1  1  te  can  undergo.  Secondly,  GPS  observations  are  not 
dependent  upon  weather  conditions  --  clouds,  high  winds,  extremes 
of  temperature,  rain  and  even  snow  do  not  significantly  hamper 
field  collection  of  data.  Host  importantly,  however,  is  the 
rapidity  with  which  GPS  data  can  be  obtained.  Whereas 
astr ogeodet 1 c  collection  may  take  several  days  (or  weeks  if  the 
weather  is  uncooperative),  a  GPS  operator  gathers  all  of  his  data, 
which  includes  a  sufficient  amount  of  redundancy,  in  only  a  few 
hours.  Furthermore,  while  it  may  take  months  to  properly  train  an 
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astrogeod.es  i  st  on  a  T-4,  a  GPS  operator  can  be  taught  his 
procedures  in  a  single  afternoon.  No  special  skills  are  required 


to  collect  GPS  data;  a  good  astrogeodes i st ,  however,  Is  a  rare 
spec  ial i st . 

A  very  noteworthy  aspect  of  GPS  data  collection  is  that  most 


of  the  work  is  performed  "in  the  lab”  and  not  in  the  field.  With 
high-speed  computers  and  highly  efficient  software,  the  more 
expensive  and  cumbersome  conventional  field  techniques  may 
eventually  become  obsolete.  As  GPS  receivers  become  more 
affordable,  many  optical  surveying  methods  will  slowly  die  out. 

In  this  paper,  I  shall  show  how  collection  of  coordinate 
difference  data  with  GPS  receivers  can  be  used  to  determine 
deflections  of  the  vertical  accurately,  efficiently  and 
1 nexpen  s i ve 1 y . 


SURVEY  GEOMETRY 


Although  there  is  an  infinite  variety  of  survey  layouts  and 
geometries,  networks  used  for  accurate  determinations  of 
deflections  of  the  vertical  will  impose  some,  although  rather 
flexible,  rules  on  the  distribution  of  the  survey  points.  One  of 
the  few  practical  applications  for  deflection  computations  is 
found  in  the  determination  of  launch  and  initial  trajectory 
characteristics  of  missiles  and  rockets.  The  inertial  guidance 
systems  of  these  space  or  sub-space  vehicles  must  be  aligned  with 
the  local  vertical  in  order  to  adequately  predict  their  paths  and 
orbital  characteristics.  Even  a  small  error  in  the  computation  of 
the  plumb  line  may  result  in  a  significant  deviation  of  the 
predicted  orbital  path  from  the  actual  one.  (Mueller,  1969; 

Kaula,  1966) 

An  optimum  survey  will  include  one  station  right  over  or  very 
near  the  point  at  which  the  deflection  of  the  vertical  is  sought. 
From  this  point  would  extend  radials  of  about  one  to  fifteen 
kilometers  in  length  to  other  peripheral  stations.  The  appearance 
of  such  a  plan  would  resemble  the  spokes  of  a  wheel.  (See  Figure 
1,  next  page. )  It  is  not  important  that  these  "spokes"  be  of 
equal  length  or  be  equally  spaced  although  an  even  distribution  is 
always  favorable. 

In  highly  benign  or  smooth  gravity  fields,  only  two 
peripheral  stations  may  be  sufficient.  As  the  complexity  of  the 
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Figure  1  . 


Plan  view  of  a  possible  survey  layout  for  determination  of 
deflections  of  the  vertical  using  GPS.  Triangle  indicates  the 
central  point  for  which  the  deflection  is  sought.  Circles  are  the 
"peripheral"  stations  to  which  GPS  measurements  are  made.  Both 
GPS  observations  and  high-order  leveling  are  required  at  all  of 
the  survey  points. 


LEAST  SQUARES  ADJUSTMENTS 


If  more  than  two  GPS  receivers  are  used  over  several  separate 
sessions  to  cover  the  survey  area,  then  the  least  squares 
adjustments  will  have  to  employ  more  complex  techniques  than  the 
ones  described  in  this  paper.  (Double  and  triple  differencing 
will  definitely  come  into  use,  for  instance.)  (Goad,  1987;  Bock 
et  al,  1986;  Beutler  et  al,  1986)  In  this  paper,  I  will  assume 
for  simplicity’s  sake  that,  in  fact,  a  radial  survey  is  conducted 
with  all  observations  on  all  survey  marks  made  simultaneously. 
Expansion  to  the  more  complex  situations  is  a  subject  already 
ue 1 1 -documented  by  many  other  authors.  (See  references) 

Using  the  notation  of  Uotila  (1986)  and  the  algorithm  of 
Sliua  (1985),  the  mathematical  model  that  ue  employ  is  one  of 
observation  equations: 

La  =  F(Xa)  =  F(X°)  +  |  X  +  ifa  I  X  + 

Xa=X°  Xa=Xl  (2) 

where  La  =  adjusted  values  of  the  observed  quantities 

Xa  =  theoretical  values  of  the  unknown  parameters. 

In  the  relative  mode  of  GPS  surveying,  one  observes 
coordinate  differences  betueen  stations:  AX(j,  AYtj  and  AZfj  ana 
eventually  solves  for  the  absolute  (though  not  necessarily 
accurate)  positions  of  those  stations:  X(-i>j,  Y  t  .  j ,  where  l 

and  j  are  the  station  identifiers.  (Sliua,  1985)  Thus,  expanding 
the  mathematical  model  from  the  generalized  form  above  yields: 


l'.VA'.li‘AV 


4 


xJa  -  X 

YJ-  - Y 

z,„  -  z 


where  station  "1"  is  the  central  point  and  all  stations  "j"  are 
the  peripheral  marks. 

One  obtains  initial  approximations  for  the  parameters  X°  by 
removing  all  known  errors  (blunders,  atmospheric  effects, 
satellite  and  receiver  clock  errors,  receiver  measurement  noise, 
etc.)  (Wei,  1985)  The  coefficient  matrix  derived  from  model 
linearization  (Sliwa,  1985)  is  then  computed  by  differentiating 
the  model  equations  with  respect  to  the  parameters: 

a  =  —  I  ( 4 ) 

°  -»X“  I 

Xa=X° 

A  vector  of  approximate  values  of  the  observables  L0  is  rigorously 


computed  by: 


L0  -  F(X°)  -  Lb 


where  Lb  is  a  vector  of  the  observed  values  of  the  observables. 
Since  the  normal  equations  can  be  represented  by: 

AJPA0X0  +  AJPL0  =0  (6) 

we  get  the  solutions  for  the  unknowns  to  be: 

XQ  -  -(A5PA0)-lATPL0  (7) 

where  P  is  a  weight  matrix.  (Uotila,  1986)  In  mu  1 1 1 -rece i ver 
simultaneous  surveying,  there  is  a  very  high  level  of  correlation 
between  observations  which  have  common  points.  Therefore, 
off-diagonal  values  in  the  weight  matrix  will  not  be  equal  to  zero 
but  in  fact  may  be  very  significant  factors  in  the  computations. 


%  *•  ■-  *.  V  \  A  A  V  .  %  .■*  S.  A  *•  -N  %  A  A 


Proper  magnitudes  for  these  correlation  coefficients  are  well 


documented  in  the  literature.  (Sliwa,  1935) 

A  set  of  adjusted  values  Xa  is  obtained  by  adding  the 
computed  solutions  for  the  unknowns  XQ to  the  initial 


approximations  X°  : 


=  x°  +  xQ 


Then  we  compute  directly: 


VIPVD  =  LjPL0  +  XjA*PL0 


where  o-o  is  the  a  posteriori  variance  of  unit  weight 

n  is  the  number  of  observations 
u  is  the  number  of  parameters 
n-u  is  the  degree  of  freedom  (Uotila,  1986) 


One  finally  computes  the  residual  vector  VD  and  the  adjusted 
values  of  the  observables  Lg  by: 

V0  -  A0Xo  +  L0  (11) 

Lg  =  F(Xg)  <  12) 

Iteration  continues  by  computing  new  A  and  L  matrices  and 
determining  better  adjusted  values,  thusly: 


A  -  ^ 


Xa=X* 


L1  =  F(X* ) 

L,  =  Ll  -  Lb 


X1  =  Xg 


and  so  forth 


The  iterative  process  of  this  adjustment  will  continue  until 


114) 


1)  X,  approaches  zero 

2)  (V(  -  Vj-,)  approaches  zero 

3  )  VfPV (  st  ab  l  1  l  zes 
4)  Lf  =  Lb  +  V,  =  F(Xf )  (Uotila,  1986) 

The  adjusted  values  of  the  observables  Lf  may  then  be  used 
in  conjuction  with  a  nominal  set  of  X,  Y,  Z  values  from  any  one 
(preferably  the  central)  station  in  the  survey  network  to  obtain 
the  X,  Y,  Z  positions  of  all  the  other  stations  in  that  network. 
It  is  very  important  to  realize  that  these  absolute  coordinates 
for  the  stations  contain  errors  far  too  great  for  our 
computations.  All  of  the  stations,  however,  have  essentially  the 
same  absolute  errors;  therefore,  their  relative  errors  are  very 
small.  We  shall  eventually  use  only  the  highly  accurate  relative 
coordinates  (coordinate  differences)  for  determination  of  the 


deflections. 


COMPUTATION  OF  THE  DEFLECTION  OF  THE  VERTICAL 


We  must  first  convert  the  absolute  cartesian  coordinates 
(X,Y,Z)  into  geodetic  coordinates  (<fr,  X,  h)  using  well  known 
transformations  presented  by  Heiskanen  and  Moritz  (1979),  Rapp 
(1984)  and  Goad  (personal  co mmun 1  cat  1  on ,  1986,1987): 

X  =  (N  +  h)  cos$  cosX 

Y  =  (N  +  h)  cos*  sinX  (  15) 

2  =  (N(l-ea)  +  h)  sin4> 

where 

N  =  7-r - =  the  length  of  the  normal  vector  from  the 

(1  -  e2  sina$)1/2 

surface  of  the  ellipsoid  to  the  minor 
axis.  (Rapp,  1984) 


Geodetic  longitude  can  be  directly  computed  by  division  of  the 
second  equation  by  the  first  (keeping  in  mind  the  quadrant  of  the 
computation  to  maintain  true  signs): 

X  =  DATAN21 Y, X )  (16) 

Computation  of  the  geodetic  latitude  <t>  and  height  above  the 
ellipsoid  h,  on  the  other  hand,,  requires  a  Newton -Rhapson 
iteration.  From  Goad  (personal  communication,  1987),  we  have: 


also 


p(4»,h)  =  p(*0,h0)  +  p$A4>  +  pxAX  +  ... 

a  P(*0)h0)  -  (N  +  h)  sin4>A4>  +  cos$Ah 


(  17) 


cos+AZ  -  siniAp 
"  N  +  h 

Ah  =  sin^AZ  +  cos^Ap 


(  18) 


where 

p  =  (N  +  h)  cost 
2  -  (N(l  -  eJ)  +  h)  3in<t> 

(  19) 

where  A4  and.  Ah  are  used  to  update  the  best  known  $  and  h 
quantities.  Initial  values  for  the  latitude  and  height  (  tQ  and 
h0  )  are  also  derived  by  Goad  (personal  communication,  1987): 


d> 

yQ 

h 


(pa  +  Z2)1/3 


(20) 


where  f  is  the  flattening  and  a  is  the  semi-major  axis  of  the  j 

reference  ellipsoid. 

t 

Once  we  have  the  absolute  geodetic  coordinates  of  all  of  the  j 

survey  points,  we  compute  the  highly  accurate  coordinate  | 

differences  in  this  same  geodetic  reference  frame.  (i.e.A4>,  AX,  Ah 
from  the  central  station  to  every  peripheral  station  in  the 

i 

network.)  As  mentioned  before,  although  the  absolute  coordinates  f 

i 

of  the  stations  themselves  are  not  acceptable,  (+  20  meters  or 
even  greater),  this  is  of  no  concern  since  we  shall  only  employ 
the  coordinate  differences  in  our  algorithm.  Even  though  the 

i 

central  station  may  have  a  large  absolute  positioning  error  of 
♦e>  he  ,  all  of  the  other  stations  measured  relative  to  this 
central  one  will  have  essentially  the  same  $B,  xo)  hB  errors.  (See 
Figure  2,  next  page.  ) 


+  *e>  X  +  XB)  h  +  he ) 


central 

station 


(♦-<  A*+4B,x+AX  +  XB,h+Ah+hB) 


A4>,  AX,  Ah 


peripheral 
station  j 


F i gur e  2. 


This  absolute  error  will  likewise  translate  into  the  N’s  but  not 
the  dN’s  of  the  survey.  I  shall  show  later  that  both  the  small 
errors  in  relative  positioning  and.  the  large  errors  in  absolute 
positioning  of  the  points  are  almost  always  negligible  in  the 
computation  of  the  deflections. 

The  geoidal  undulation  for  each  point  can  be  computed  using 
the  difference  between  the  geometric  height  above  the  ellipsoid 
(h)  and  the  orthometric  height  above  the  geoid  (H): 

N,  =  h,  -  H,  (21  ) 

where  the  or thouetr ic  heights  can  be  obtained  by  conventional 
high-order  leveling  techniques. 

Because  we  only  have  highly  accurate  coordinate  differences 
between  the  stations,  we  must  not  use  the  absolute  geoidal 


At  this  point,  ue  have  very  good  &♦,  AX,  Ah  and  dNt  j’s  for  aH 


of  the  lines  in  the  network.  Given  only  this  data,  ue  can  model 
the  surface  of  the  geoidal  differences  using  interpolation  by 
least  squares  fitting  of  polynomial  surfaces.  The  specific 
algorithms  that  ue  employ  will  be  dependent  upon  the  number  of 
survey  points  in  our  netuork.  According  to  Lancaster  (1986),  the 
general  form  of  a  quadratic  function  in  tuo  variables  i  s  •• 

p(x,y)  =  a,  +  a2x  +  a3y  +  a4x2  +  a5xy  +  a6y2  (24) 

If  f(  is  the  geoidal  undulation  difference  between  stations,  then 
we  want  to  adjust  the  function  p  by  choosing  the  constants  to  a6 
so  as  to  minimize 

E(p)  =  £  (p(x,,y,)  -  fj)2  (25) 

1  =0 

(Lancaster,  1986) 

This  function  E  will  have  a  minimum  only  when  -  =  0  for  i  =  1  to 

*a, 

6  and  will  yield  six  linear  equations  as  functions  of  the  a’s. 

The  six  normal  equations  are: 

(/V  +  l)d[  +  (2x,)a2  +  (2y,)fl3  +  (2xj"  )a4  (26) 

+  (Zxi>'i)a5  +  (2y2K  =  2/,-, 

(Z.v,)a,  +  (2x2)a:  +  (2x,y,)a3  +  (2x3)a4 
+  (. 2x2y(>5  +  (2x,y2K  =  1.x,  f, , 

(SyJaj  +  {lx,y,)a2  +  (I  yj)a3  +  (2x2y,>4 
.  +  {lx,yj)a}  +  (2y3)a6  =  2y,/,, 

(2  xj)a{  +  (2  x])a2  +  (2x2y,)a3  +  (2x-)a4 
+  (2  x]y,)a5  +  (2  xjy;)a6  =  2  x2/,, 

(2x,y,)al  +  (2x2y>2  +  (2  x,y;)a3  +  (2x3v,)a4 
+  (2  xfyf)a5  +  (2x,y3  )a6  =  2  x,y,/,, 

(2y2)a,  +  (2x,yr)a:  +  (2y?)a3  +  (2  x;y;)aA 
+  (2x,y3)a5  +  (2  yf)a6  =  2  y2/,, 


and  the  vectors  : 


Si 

'  AN(40,\a)  ‘ 

£N 

a2 

*N 

a3 

a\ 

— 

a2N 

a« 

2 

a2N 

a5 

a2  N 

a6 

L  J 

dN, 

dNa 

.  dNn 


(  30  ) 


The  normal  equations  then  take  on  the  form: 

VTV3  =  V T? 


(31  ) 


or,  to  facilitate  a  solution: 

t  =  (VTV)~‘VT?  (32) 

As  a  warning,  one  should  be  very  aware  of  the  possibility  of 

creating  i 1 1 ~cond i t i oned  VTV  matrices  caused  by  factorization  of 
the  coefficient  matrix  or  by  a  c 1 ose-to-orthogona 1  (pathological) 
survey  network.  (Lancaster,  1986)  If  the  network  must  be  of  an 
orthogonal  design,  then  one  would  advance  to  a  more  complex 
interpolation  routine  involving  cubic  or  bicubic  splines. 
(Lancaster,  1986) 

An  advantage  of  the  above  system  is  that  one  may,  if  one  has 
sufficiently  accurate  information,  impose  a  weight  matrix  (P)  on 
the  normal  equations  to  obtain: 

t  ~  (VT PV)-‘VTP?  (33) 

The  diagonal  coefficients  in  this  weight  matrix  could  be  dependent 
upon  the  distances  between  points,  the  accuracy  with  which  various 


& 


parts  of  the  field  surveys  are  conducted,  recency  of  the  collected 


data,  and  so  forth. 


In  this  paper,  however,  I  shall  assume  no 


such  knowledge  and  constrain  the  ?  matrix  to  the  identity  matrix. 


Interpolation  with  three  stations: 


The  number  of  survey  points  observed  will  determine  which 
case  of  the  polynomial  should  be  used.  If  we  only  observe  two 
peripheral  stations  (for  a  total  of  three  survey  points),  then  we 
cannot  use  the  polynomial  form  discussed  above  but  must  rather 
revert  to  a  highly  shortened  version  of  it.  With  three 
non-co 1 1 i near  survey  points,  we  can  only  construct  a  plane  to 
approximate  the  undulation  difference  surface.  (See  Figure  3.) 
The  equation  of  such  a  plane  is 

dN(A4>,  AA)  =  AN(4>0,A0)  +  ^  A*  +  ~  AX  (34) 

where  4>0IA0  and  AN($01Ao)  are  undefined.  (Thus  we  can  set  them  a 
priori  equal  to  zero.)  According  to  Lancaster  (1986),  we  first 
compute  A  where 


A  =  det 


♦r 


1 

1  A*, 

1  A*2 


^  O 

AA, 

AA2 


=  A<t> ,  AAa  —  A*2AA, 


then  we  get  directly 


<*N 


=  A* 


det 


!*(♦„, A„) 
dN, 
dNa 


AA, 
AA  2 


and 


3- »- 


det 


A*, 


N(*0,A0) 

dN, 


=  A~ 1  [dNa*,  -  dN, *2 ] 


(  35  ) 


=  A~ 1  [dN, AAa  -  dNa AA, ]  (  36  } 


(  37  ) 


Figure  3. 


Plan  view  of  a  survey  network  containing  only  two  peripheral 
stations  (circles).  In  this  case,  the  geoidal  undulation 
difference  surface  can  only  be  approximated  by  a  plane. 


<£> 

Figure  4 , 


Survey  network  containing  three  peripheral  stations  (for  a  total 
of  four  survey  marks).  The  geoidal  undulation  difference  surface 
will  be  approximated  by  a  curved  surface  derived  by  least  squares 
adjustments  of  polynomials.  Similar  layouts  will  be  made  for  four 
or  more  peripheral  stations. 


1  v' 


Interpolation  with  four  stations: 


For  a  situation  involving  three  peripheral  points  (four 
stations  in  all  —  see  Figure  4.),  ue  cannot  (usually)  approximate 
the  surface  as  a  plane.  We  can,  however,  use  a  more  advanced 
version  than  the  one  above.  Once  again  ue  have: 

dN(A<fr,AX)  =  AN(*0,X0)  +  T7  i*  +  AX  (  38  ) 

but  in  this  instance,  ANOo ,  XD)  cannot  be  set  a  priori  equal  to 
zero.  The  V  matrix  that  we  described  before  will  be  a  3x3  of  the 


form 


V 


1  A$i  AX! 

1  A$a  AX  2 

1  a*3  ax3 


and 


3  - 


AN(*0,X0)  ' 
£N 

d  ($> 


aN 

«x 


dN, 

dN2 

dN3 


Here,  ue  can  immediately  compute 


Once  again. 


t  =  (VTV)-lVT? 

ii]  s  W 

—  and  — -  will  be  directly  resolved. 


(  39  ) 


(  40  ) 


(4  1) 


Interpolation  with  five  or  more  stations: 


Finally,  there  is  the  most  generalized  case  in  which  we 


observe  four  or  more  peripheral  stations  in  addition  to  the  one 


central  station.  Now,  the  equation  of  the  surface  to  be  modeled 


is,  as  before: 


«2N 


« 2N 


dN(  A4> ,  AX )  =  AN(*0,X0)  +  —  At  +  —  AX  +  ^  A^2  +  A*AX  +  —  AX: 


(  42  ) 


and  the  matrices  V,  f  and  a  are  of  the  same  form  as  in  equations 
29  and  30.  The  solution  to  this  model  is: 

3  =  (VTV)"*VT?  (  43) 


4  N  8  ^ 

where  —  and  —  can  be  directly  solved.  The  additional  terms 

•  ▼  *  A 


<»JN 
a  4>2  ’  a$a\ 


2N  a  2N 

and  will  not  be  necessary  in  the  finai  computation 


would  otherwise  be  erroneously  dumped  into  —  and  — . 

*♦  *X 


of  the  deflections,  but  they  will  "absorb"  minor  deviations  that 

(This 

argument  may,  of  course,  be  taken  to  extremes  by  computing  still 
higher  derivatives  of  the  geoidal  undulations  if  even  more  points 


are  observed,  but  these  corrections  will  truly  be  inf mitesmal . ) 

a  N  a  fj 

Once  —  —  and  an  average  latitude  $  for  the  survey 

<*9  a\ 

network  are  all  computed,  the  determination  of  the  deflection  of 


the  vertical  components  is  very  simple.  From  Heiskanen  and  Moritz 

(1979),  we  have  in  a  north-south  direction: 

.  _  dN  1  <»N 

*  "  e*  dS$  K.  a +  (44) 

in  an  east-west  direction: 

dN  _  1  «N 

v  *X  dSx  Rcos4>  a\  (45) 

One  may  use  algebra  and  trigonometry  to  obtain  the  total 


deflection  magnitude  and  azimuth: 

£  =  7  ^ 2  +  v2  (magnitude)  (46) 

<xt  =  DATAN2(  £  ,  tj)  (azimuth)  (47) 


„■>'  .">'  '."*■  ^:,^^^>^jr’jr<in'±-*jr’y-'y~*j,-’j.- ~ ->  *jrw)r-jr-'jf.  nxn/nA’V  .'vcwv» 
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ERROR  ANALYSIS 

In  order  to  confirm  the  validity  of  the  arguments  made  here 
and  to  check  the  accuracy  of  the  deflections,  a  thorough  error 
analysis  must  be  performed  on  results  obtained  using  this  GPS 
technique.  Unfortunately,  a  last-mi  uite  failure  in  one  of  the 
Geodetic  Science  Department’s  Kaypro  2000  computers  (which  is  used 
in  conjunction  with  the  Trimble  model  GPS  receiver),  prevented  any 
collection  of  field  data.  For  this  reason,  GPS  and  leveling  data 
collected  and  adjusted  for  another  project  must  be  used  in  our 
analysis.  Although  not  completely  valid,  these  tests  should  at 
least  provide  some  idea  as  to  the  accuracies  one  might  expect  if 
one  were  to  obtain  actual  GPS  data  for  the  determination  of  the 
deflection  of  the  vertical. 

The  consensus  in  current  literature  lists  the  accuracy  of 
GPS-determined  coordinate  differences  to  range  from  1  to  3  ppm. 

(Goad,  personal  communication,  1S87)  Goad  (1987)  has  used  a 
formula  of  3mm  +  lxlO'6  x(L)  for  computing  the  error  in  baselines 
of  length  L.  The  3mm  factor  comes  from  "noise"  in  the  system.  An 
"optimum"  length  baseline,  however,  is  sti 1  *  a  debatable  concept 
and  should  not,  at  least  for  now,  be  specified.  As  mentioned  in 
the  survey  layout  section  of  this  report,  baselines  from  1  to  15 
kilometers  should  be  sufficient  to  adequately  model  the  geoidal 
undulation  difference  surface.  (Lines  of  1  to  5  kilometers  in 
length,  however,  are  probably  most  appropriate  for  our 
interpolation  schemes. ' 


A  value  must  be  placed  on  thp  desired  accuracy  for  the 


deflection  of  the  vertical.  Current  astrogeodet i c  techniques, 
using  Wild  T-4  astrotheodo 1 i tes ,  are  capable  of  reaching  values  in 
deflections  of  ±  .5  seconds  of  arc.  (Mueller,  personal 
communication,  1987)  GPS  methods  must  therefore  accomplish  at 
least  this  accuracy  as  well.  (We  recall  that  the  purpose  of  using 
GPS  for  determinations  of  deflections  of  the  vertical  is  not  to 
increase  accuracy,  but  rather  to  increase  efficiency  and 
cost-effectiveness  of  making  such  computations,  therefore  ue  are 
not  searching  for  better  results.) 

Using  GPS  data  collected  in  the  Columbus,  Ohio  area  by  the 
Franklin  County  Engineers  Office  in  March  1987  and  adjusted  / 
analyzed  by  Dr.  Goad’s  multi-receiver  TRIMVEC  program  (WGS-84 
parameters),  ue  have: 


Name 


Elevation 


Station  A  (Central  station)  Britton 

Station  B  (Peripheral  station)  18-83 

Station  C  (Peripheral  station)  Clark 


f  .  59  m 
2  -  l .  7  4  m 
273 . 50  m 


For  line  A-B  (Britton  to  18-83),  we  compute  the  ueighted  mean 
(from  two  observation  sessions): 


d*  =  -0.089487058  deg 

dX  =  0.094818595  deg 

dh  =  -54.0591445  m 
dH  =  -53.85  m 

dN  =  -  0.2091  m 


.* 

I] 

1*1 


,•41 


f 

f 


F 

I 

R 


a 


For  line  A-C  (Britton  to  Clark),  ue  compute  (from  one  observation 
session) : 

=  -0.000158510  deg 

dA  =  0. 158519430  deg 

dh  =  1 . 6263  m 

dH  =  1. 910  m 

dN  =  -0.2837  m 


Since  ue  have  only  three  survey  points,  ue  must  approximate  the 
geoidal  undulation  difference  surface  with  a  plane.  Using 
Lancaster’s  methods: 


A  =  det 


1  0 
1  -.089487058 

1  -.000158510 


0 

.  094818595 
.  158519430 


-.  014170408 


A-1  =  -70.56959962 


and 

1  0 
1  -.2091 

1  -.2837 

and 


£N 

a  4,  ~  A 


det 


0 

.  094818595 
.  158519430 


£N 

*X 


1  0 


A  1  det 


-.089487058 


1  -.000158510 


0 


-.2091 


-.2837 


.  440804353 


-1. 789245193 


computing  directly,  ue  obtain: 


V 

total  magnitude 
azimuth 


-0. 82  seconds 
4 . 33  seconds 
4.41  seconds 
100°  43’  24‘L  5 


Rapp  (1987,  personal  communication),  has  obtained  the  following 
results  using  a  360  degree  model  of  the  geoid  (minimum  wavelength 
of  50  km ) : 


Station 


£  seconds 


v  seconds 


A 

B 

C 


1. 36 
1. 26 
2. 03 


2.  57 
1 . 90 
2.  18 


average 


1 . 55 


2.22 


Interpolation  off  a  map  of  the  geoid  in  Central  Ohio  (Engelis, 

Rapp)  (See  Figure  5),  yeilds: 

£  =  -0.5  seconds 

ij  =  4.0  seconds 

Although  there  is  little  correlation  between  GPS  derived  values 
and  those  obtained  with  Rapp’s  gravity  model,  (due  mostly  to  the 
inability  of  Rapp’s  model  to  resolve  deflections  from  wavelengths 
of  less  than  50  km  in  length),  there  is  an  extremely  high  correla¬ 
tion  with  the  values  obtained  by  interpolation  from  the  geoidal 
m-j.  The  deflections  computed  with  this  GPS  technique  therefore 
appear  to  be  valid  and  reliable  (at  least  for  this  particular  test 
area.  ) 

The  next  step  is  to  subject  our  GPS  method  to  various  errors 
(both  absolute  and  relative)  and  see  how  they  propagate  into  our 


solut ion . 


Error  gjj 

Impose  an  absolute  positioning  error  on  Station  A  (Britton) 
of  +20  meters  in  the  X  direction.  (Stations  B  and  C  will  of 
course  experience  the  same  shift.): 


£N 

a  4,  ~ 

. 415969848 

m 

£N  _ 
a\ 

-  1 . 5217939  10 

m 

(  = 

-0.  77 

sec 

V  = 

3.  68 

sec 

Error  #2: 

Impose  a  relative  positioning  error  on  stations  B  and  C  of 


(  lx  10' 6 )  x  ( length 

0  f 

line 

)  : 

Errors 

on 

line 

A-B: 

x : 

.  009 

m 

y : 

-.005 

m 

z  : 

.  008 

m 

Errors 

on 

line 

A-C  : 

x  : 

.013 

m 

y : 

-.002 

m 

z  : 

.  000 

m 

£N 

. 

495630234  m 

<»N  _ 

-1 . 

724844248  m 

i  = 

-0. 

92 

sec 

v 


4.  17 


sec 


I  . 

♦i 


r*ii ijy  *ii  *it  *d 


W/jUASVt 

■jl 


Error  tt3: 

Impose  an  error  on  the  absolute  elevation  above  the  geoid 
of  point  A.  First  order  standards  require  that  elevations  be 
measured  to  a  tolerance  of  3mm  x  (distance):  where  the  distance 

is  in  kilometers.  (FGCC,  1980)  The  maximum  error  ue  may 
experience  is  then: 


1 ine  A-B  AH  =  .011  m 


line  A-C  AH  =  .011  m 


So  ue  get  : 


.246832384  m 


=  -1.856307937  m 


-0.46 


4. 49 


In  tabular  form:  (all  quantities  are  in  seconds) 


No’  Errors 


Abs.  Error 


Re  1 .  Error 


Lev.  Error 


We  see  that  in  almost  every  case,  the  errors  ue  imposed  did 
not  deviate  significantly  from  the  original  values.  If  our 
tolerance  is  +.5  arcseconds  then  all  but  the  absolute  error  in  the 
meridian  direction  are  thoroughly  acceptable.  Overall,  the 
results  are  characterized  by  a  striking  consistency. 


CONCLUSIONS 


I 


The  similarity  between  the  results  obtained  with  the  GPS 
algorithm  described  in  this  paper  and  those  from  Dr.  Rapp’s  map  of 
the  geoid  imply  that  the  use  of  GPS  for  determining  deflections  of 
the  vertical  is  a  valid  procedure.  Of  course,  more  data  and  more 
comparisons  are  required  for  even  stronger  confirmation.  There  is 
no  theoretical  reason  why  the  GPS  methods  could  not  work 
exceptionally  well  in  benign  geoidal  terrains.  In  regions  of  high 
variability,  however,  one  must  be  very  careful  with  one’s 
conclusions . 

The  test  bed  that  I  used  in  my  analysis  was  locally  rather 
benign  but  regionally  quite  variable.  Since  the  isarythms  of  the 
geoidal  surface  were  only  slightly  curved  in  the  test  area,  we  can 
be  reasonably  sure  of  the  quality  of  the  results.  Had  the 
variability  been  more  extreme,  either  more  points  or  shorter  lines 
would  have  been  required.  It  seems  as  though  in  the  future  one 
will  need  at  least  some  idea  of  the  characteristics  of  the  geoidal 
surface  prior  to  planning  a  survey  layout.  With  modern  satellite 
methods  and  worldwide  gravity  models,  however,  this  should  not  be 
a  major  stumbling  block. 

Overall,  I  would  use  GPS  methods  to  determine  deflections  of 
the  vertical  in  areas  which:  1)  have  benign  geoidal  terrains;  2) 

commonly  experience  poor  weather;  3)  cannot  be  occupied  for  more 
than  several  hours  (including  hostile  and  unsafe  regions);  4) 


require  immediate  data;  and  5)  must  be  surveyed  as  cheaply  as 
possible.  I  would  not  use  GPS  methods  in  areas  which:  1)  have 
extremely  variable  geoidal  undulations;  and/or  2)  experience  a 
high  amount  of  radio  frequency  interference  (RFI).  Overall,  the 
positive  aspects  of  GPS  surveying  seem  to  outweigh  the  negative 
ones.  Determinations  of  deflections  of  the  vertical  using  GPS 
will  probably  become  the  method  of  choice  for  many  agencies  in  the 


near  future. 
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The  following  pages  are  output  from  Dr.  Goad’s  mu  1 1 i rece  i  ver 
Triivec  program  which  adjusts  GPS  data.  Three  sessions  are 
i nc 1 uded : 


1 . 

Britton 

-  Clark 

5 

March 

1987 

2. 

Britton 

-  18-83 

Session 

A 

6 

March 

1987 

3. 

Britton 

-  18-83 

Session 

B 

6 

March 

1987 

Default  atmospheric  values  and  adjusted  antenna  heights  were 
used  in  these  programs.  All  data  was  collected  by  the  Franklin 
County  Engineers  Office.  (Point  of  contact  is  Mr.  Dean  Ringle.) 
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TRIMBLE  NAVIGATION,  LTD 
585  NORTH  MARY  AVENUE 
SUNNYVALE,  CALIFORNIA  94086 

U'S'A'  8  il  ITFC'V  —  Ll  KK 

PROGRAM  TRIMVEC 

GPS  RELATIVE  POSITIONING  SOLUTION 
VERSION  87. 106MB 

Filename:  a : br i tc lar . a64 
Coordinate  system  -  WGS-84 

Type  solution:  Double  difference 

Start  date/time:  1987/  3/  5  6:58:30.  day  of  year 

Stop  date/time:  1987/  3/  5  8:26:15.  day  of  year 


64  tow  370710 
64  tow  375975 


Data  avai lable : 

station:  1 

sat :  3  ! . 

sat :  9  1 . 

sat :  1 1  I . 

sat :  6  1 . 

sat:  13  ! 


station:  2 

sat :  3  !  .  .  . 
sat :  9  !  .  .  . 
sat :  1  1  I . . . 
sat:  6  i... 
sat : 13  ! 


Ephemens  file  used:  a  :  br  i  ta064  .  eph 


SATELLITE 

AODE ( hr .  ) 

HEALTH  WEEK 

NO. 

TOW (sec)  URA(m) 

3 

2.84 

0 

373 

370590.00 

8.0 

9 

6.83 

0 

373 

370590.00 

8.0 

1  1 

3.98 

0 

373 

370590. 00 

8.0 

6 

4.55 

0 

373 

370620.00 

4.0 

13 

8.53 

0 

373 

372660.00 

8.  0 

Broadcast 

sate  1 1 i te  c lock 

correct  ion 

values 

prn 

afO 

af  1 

af  2 

toe 

3  - 

. 3044493496D-05 

568434 1886D- 12 

. 0000000000D+00 

. 3744D+06 

9 

. 3178184852D-04 

.  165982783  ID- 10 

- . 2775557562D- 16 

. 3744D+06 

1  1 

. 603 18969 19D-04 

51 15907697D- 1 1 

. 0000000000D+00 

. 3744D+06 

6  - 

. 34043379 13D-03 

-.  1546140993D-10 

.  0000000000D+00 

. 3744D+06 

13 

. 9765336290D-04 

. 22737 36754 D- 1  1 

.  0000000000D+00 

. 3780D+06 

Message  file  for  station  1 


Station  ID:  BRITA064 
Date:  MAR.  4,1987 


Sess ion :  064A 

Location:  BRITTON  RD.  HILLIARD 

Observer’s  name:  BRIG 

Receiver  serial  #:  4871 

Antenna  serial  #:  0117 

External  computer  serial  #:  864C 

Antenna  height:  1.3770  meters  (54.213  inches) 

BRITTON  RD.  HILLIARD 

STA.  BRIT 

SESSION  A 


■  tvtvvvvvv*  vnr">'’_>  ->  ■ 


-  uncorrected  to  the  vertical 


Origin  of  station  1  coordinates  :  Best  C/A  code  tracking  solution 
STATION  (mark)  1 


input  data  file  1 
antenna  height(m) 

met  values  used:  i 


L  :  a : br ita064. dat 
1  1 . 343 

pressure(mb)  1010.0 
temperature (deg  C)  20.0 
relative  humidity(%)  50.0 


X 

(m) 

584454.502 

lat 

( dms  ) 

N 

40 

3 

10. 19823 

y 

( m  ) 

-4854000.441 

e  Ion 

( dms ) 

E 

276 

51 

56.69108 

z 

(m) 

4082582.953 

ulon 

( dms  ) 

W 

83 

8 

3. 30892 

ht  ( m ) 

163.5942 

Message  file  for  station  2 


Observer’s  name:  KLS 

Receiver  serial  #:  4874 

Antenna  serial  it :  0116 

External  computer  serial  #:  5501905C 

Antenna  height:  1.5926  meters  (62.701 

CLOUDY  AND  COLD.  THE  STATION  OCCUPIED 

BE  USED  AS  PHOTO  CONTROL. 


inches)  -  uncorrected  to  the  vertical 
IS  STATION  CLARK  WHICH  WILL 


STATION  (mark)  2 

input  data  file  1 
antenna  height(m) 


a : c 1 ara064 . dat 
1. 563 


met  values  used: 


pressure(mb)  1010.0 
temperature ( deg  C)  20.0 
relative  humidity(%)  50.0 


X 

( m ) 

597883.283 

lat 

( dms  ) 

N 

40 

3 

9.62758 

y 

( m ) 

-4852377.340 

e  1  on 

( dms  ) 

E 

277 

1 

27.36113 

z 

( m ) 

4082570. 527 

ulon 

( dms  ) 

W 

82 

58 

32. 63887 

ht  ( m  ) 

165.2205 

Station  ID:  CLARA064 
Date:  5  MARCH  1987 
Session:  064A 

Location:  DESALLES  HIGH  SCHOOL  KARL  ROAD  FRANKLIN  COUNTY  STATION  CLARK 


dz (01  ) 


Vector  covariance  matrix  (cm**2) 
dx(01)  dy(01) 


dx(01)  .0030774 

dy(01)  .0032010  .0203249 

dz(01)  -.0022806  -.0124487 

Vector  correlation  matrix  : 


. 0090608 


dx  ( 0 1 ) 

dy  (01  ) 

dz (01  ) 

dx(01)  1. 

OOOOOOO 

djr(Ol) 

4047446 

1. OOOOOOO 

dz(01)  -. 

4318892 

-.9173319 

1.0000000 

STATION  1  TO 

STATION 

2 

Rdop  is 

.  107680 

slope  distance 

( m ) 

13526.5208  sigma 

( cm ) 

064 

normal  section 

az imuth 

(dms)  90 

1  24.78 

vertical  angle 

( dms  ) 

0 

-3  13.61 

east(m)  north(m) 

up  ( m ) 

13526.514 

-5. 

560 

Delta  lat(dms) 

0 

0  -.57064 

Delta  lon(dms) 

0 

9  30.67005 

Delta  ht ( m ) 

1 . 6263 

corre lat ions 

for  baseline  1: 

dx 

dy 

dz 

trop 

bias  1 

bias  3 

bias  4 

dx  1 . 

OOOOOOO 

dy 

4047446 

1. OOOOOOO 

dz 

4318892 

-.9173319 

1 . OOOOOOO 

trop 

OOOOOOO 

. OOOOOOO 

. OOOOOOO  ] 

l.  OOOOOOO 

bias  1 

OOOOOOO 

. OOOOOOO 

. OOOOOOO 

. OOOOOOO 

1 . OOOOOOO 

bias  2 

OOOOOOO 

. OOOOOOO 

. OOOOOOO 

. OOOOOOO 

. OOOOOOO 

bias  3 

OOOOOOO 

. OOOOOOO 

. OOOOOOO 

.OOOOOOO 

. OOOOOOO 

1. 

OOOOOOO 

bias  4 

OOOOOOO 

. OOOOOOO 

. OOOOOOO 

. OOOOOOO 

. OOOOOOO 

OOOOOOO 

1 . OOOOOOO 

Solution 

S  i  gma 

dx  ( m ) 

13428.781 

.  001 

dy  ( m ) 

1623.101 

.  001 

dz  ( m ) 

-12.426 

.  001 

trop  (%) 

.000 

.  000 

bias  1  (cycle ) 

.000 

.  000 

bias  2  ( cycle ) 

.000 

.  000 

bias  3  (cycle) 

.000 

.  000 

bias  4  (cycle) 

-14.000 

.  000 

-12.697 


bias  2 


1.  OOOOOOO 
.  OOOOOOO 

.  OOOOOOO 


Results  of  integer  bias  search: 

.81201  15.55343 


28.66430 


0 


>  i.vav'iviuvi'Ai.'iiAbU 


AsmuiVn'iVj 


prn  9  0 
prn  i  1  0 
prn  6  0 
prn  13  0 

Ratio  sui-o f -square s ( 2 )  to  suo-of  squares!  1) 


Interval  between  epochs  (sec)  15 
Epoch  increment  1 

Number  of  measurements  used  in  solution 
Number  of  measurements  rejected 
RMS  (cycles)  .033 

Elevation  mask  (deg)  15.0 

Edit  multiplier  3.5 

Modified  Hopfield  troposphere  model  used 


19.  15 


4 


A 

5 


Best  tracking  C/A  code  positions 
Station  1 
Pdop  15.  1 


x  (m) 

584454.625 

1  at 

( dms ) 

N 

40  3 

10. 19823 

y  ( m ) 

-4854001. 462 

e  Ion 

( dms ) 

E 

276  51 

56.69108 

z  ( m ) 

4082583.817 

ulon 

( dms ) 

W 

83  8 

3.30892 

ht  ( m ) 

164.9372 

c lock  of  f  set ( s  ) 

. 67806059D 

-03 

freq  off  set ( s/s ) 

.  16051 102D 

-10 

Code  cal ibration (m) 

Carrier  ca 1 i brat i on ( m ) 

1  -  2 

-. 5271 

-. 

0009 

1  -  3 

-. 1090 

. 

0005 

1  -  4 

-.3663 

-. 

0006 

1  -  5 

-1.0829 

0016 

Stat ion 

2 

Pdop 

12.7 

x  ( m ) 

597904.922 

1  a  t 

( dms  ) 

N 

40  3 

9.51469 

y  ( m ) 

-4852390.070 

e  1  on 

( dms  ) 

E 

277  1 

28.20151 

z  ( m ) 

4082578. 823 

ulon 

( dms ) 

W 

82  58 

31. 79849 

c  lock  off  set ( s ) 
freq  of  f  set ( s/s ) 


Code  ca 1 i brat  ion ( m ) 
1-2  -.1251 
1  -  3  2.2444 
1  -  4  -.5557 
1  -  5  2.  1747 


ht  (m)  182.2562 

-. 77472944D-01 
. 60433 1 10D-08 

Carrier  ca 1 i brat  ion ( m ) 

.  0007 
.  0005 
.  0005 
.  0007 
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TRIMBLE  NAVIGATION,  LTD 
585  NORTH  MARY  AVENUE 
SUNNYVALE,  CALIFORNIA  94086 
U.  S. A. 

PROGRAM  TRIMVEC 

GPS  RELATIVE  POSITIONING  SOLUTION 
VERSION  87. 106MB 


tifliTTCA,1  - 

a/' 


iS'tS 

A 


File'name:  a : br i t 1883 . a65 
Coordinate  system  -  WGS-84 

Type  solution:  Double  difference 

Start  date/tine:  1987/  3/  6  ( 

Stop  date/time:  1987/  3/  6  i 


Data  available: 


stat ion 
sat :  3 


6:52:30.  day  of  year  65  tow  456750. 

8:19:45.  day  of  year  65  tou  461985. 


sat :  11 
sat :  13 


stat ion 


Ephemeris  file  used:  a : br i t a065 . eph 


SATELLITE 
3 
6 
9 
1  1 
13 


prn 

3 
6 
9 
1 1 
13 


DE  (  hr  .  ) 

HEALTH  WEEK 

NO. 

TOW ( sec  ) 

URA ( m ) 

2.84 

0 

373 

456630.00 

5.  7 

5.69 

0 

373 

456660. 00 

4.0 

6.83 

0 

373 

456630. 00 

8.  0 

10.24 

0 

373 

456630.00 

8.0 

8.53 

0 

373 

458820.00 

5.  7 

ite  clock 

correct i on 

values 

afO 

3090 128303D-05 
34  17823464D-03 
33 16905349D-04 
5984678864D-04 
978438 1837D-04 


af  1 

-.5684341886D-12 
.  15916 15728D-  10 
. 16029844 12D-10 
. 545696821  ID- 1 1 
. 2 1600499  17D-  1  1 


af  2 

0000000000D+00 
OOOOOOOOOOD+OO 
2775557562D- 16 
2775557562D- 16 
0000000000D+00 


toe 

. 4608D+06 
. 4608D+06 
. 4608D+06 
. 4608D+06 
. 4644D+06 


Message  file  for  station  1 


Station  ID:  BRITA065 
Date :  MAS. 6,  1987 


Location:  BRITTON  SB.  HILLIARD 
Observer’s  name:  BRIG 
Receiver  serial  #:  4871 
Antenna  serial  #:  0117 
External  computer  serial  #:  864C 

Antenna  height:  1.3560  meters  (53.386  inches)  -  uncorrected  to  the  vertical 
BRITTON  RD.  HILLIARD 
STA.  BRIT 
SESSION  A 


Origin  of  station  1  coordinates 
STATION  (mark)  1 


Best  C/A  code  tracking  solution 


input  data  file  1 
antenna  height(m) 


a : br ita065. dat 
1 . 321 


met 

values 

used:  pressure(mb) 

10 

>10.0 

temperature ( deg  C 

) 

20.0 

re lat i ve 

humidity 

(%) 

50.0 

X 

( m ) 

584396.324 

1  at 

( dms ) 

N 

40 

3 

10.64840 

y 

( m ) 

-4854051. 276 

e  Ion 

( dms ) 

E 

276 

51 

53.99783 

z 

( m ) 

4082637.675 

wlon 

( dms  ) 

W 

83 

8 

6.00217 

ht  ( m 

) 

232. 

1 165 

Message  file  for  station  2 


Station  ID:  I883a065 
Date:  03-06-87 
Session:  065 

Location:  mon  1883  in  park  at  corner  of  central  6<  mckinely 

Observer’s  name:  dip 

Receiver  serial  #:  4872 

Antenna  serial  # :  2704a00128 

External  computer  serial  #:  5501536c 

Antenna  height:  1.6080  meters  (63.307  inches)  -  uncorrected  to  the  vertical 
warm  40  degrees  roughly  &  very  calm 


STATION  (mark)  2 


input  data  file 
antenna  height 

met  values  used 


1  : 

a : 

1883a065. dat 

) 

1.579 

pr 

ess 

ure ( mb ) 

1010. 

0 

te 

mpe 

rature (deg  C  ) 

20. 

0 

re 

lat 

ive  humidity(%) 

50. 

0 

1  at 

( dms  ) 

N 

39 

57 

48.49478 

e  Ion 

( dms ) 

E 

276 

57 

35.34498 

wlon 

( dms  ) 

W 

83 

2 

24.65502 

ht  ( m 

) 

178.0629 

x  (m) 
y  ( m ) 
z  ( m ) 


593197.534 

-4859377.496 

4074992.037 


l  X  VX  '-V  i  «  \ 


N.  *«  ■  V.  \  *v  •  <-  »C 


Vector  1  originates  at  station  1  ends  at  station  2 


Vector  covariance  matrix  (cm**2) 
dx(01)  dy(01) 


dz ( 01  ) 


dx(01)  .0019437 

dy(01)  .0024505  .0075512 

dz(01)  -.0020715  -.0047747  .0042756 

Vector  correlation  matrix  : 

dx(Ol)  dy(01)  dz(01) 


dx(01)  1.0000000 

dy(01)  .6396375  1.0000000 

d± (01)  -.7185848  -.8403219  1.0000000 


STATION  1  TO  STATION  2 


Rdop  is  .095704 

slope  distance  (m)  12817.3986  sigma  (cm)  .041 


normal  section  azimuth  (dms)  140  47  47.44 
vertical  angle  (dms)  0-17  57.31 


east(m)  north(m)  up ( m ) 


8101.466 


-9932.  144 


-66.944 


Delta  1  at ( dms )  0  -5  22.15362 
Delta  lon(dms)  0  5  41.34716 
Delta  ht(m)  -54.0536 


correlations  for  baseline  1 


dx 

dy 

dz 

trop 

bias  1 

bias  2 

bias  3 

bias  4 

dx 

1. 0000000 

dy 

. 6396375 

1. 0000000 

dz 

-.7185848 

-.8403219 

1 . 0000000 

trop 

. 0000000 

. 0000000 

.  0000000 

1 . 0000000 

bias  1 

. 0000000 

. 0000000 

.  0000000 

. 0000000 

1 . 0000000 

bias  2 

. 0000000 

. 0000000 

.  0000000 

.  0000000 

. 0000000 

1. 0000000 

bias  3 

. 0000000 

. 0000000 

.  0000000 

. 0000000 

. 0000000 

.  0000000 

1. 0000000 

bias  4 

. 0000000 

. 0000000 

.  0000000 

. 0000000 

. 0000000 

. 0000000 

. 0000000 

1. 0000000 

Solution 

Sigma 

dx  ( m ) 

8801. 21 1 

.  000 

dy  ( m ) 

-5326.220 

.  001 

dz  ( m ) 

-7645.638 

.  001 

trop  (%) 

.  000 

.  000 

bias  1  (cycle) 

.000 

.  000 

bias  2  (cycle) 

.  000 

.  000 

bias  3  (eye  le  ) 

.000 

.  000 

bias  4  (cycle) 

.000 

.  000 

Results  of 

integer  bias  search: 

4.07169 

16 

.4  1816 

25. 

36119 

prn  6 


>  /i  .%  .% 


TO* 


TO 


TOC 


-  P'r 

prn  11  0  0 
prn  13  0  0 
Ratio  sum-of -square s ( 2 )  to  sum-of-squares (  1  ) 


Interval  between  epochs  (sec)  15 
Epoch  increment  1 

Number  of  measurements  used  in  solution 
Number  of  measurements  rejected 
RMS  (cycles)  .025 


4.03 


1012 

172 


Elevation  mask  (deg)  15.0 

Edit  multiplier  3.5 

Modified  Hopfield  troposphere  model  used 

Best  tracking  C/A  code  positions 
Station  1 


1 

■ 

Pdop 
x  ( m ) 

14.3 

584396. 445 

1  at  ( dms  ) 

N 

40 

3 

10. 64840 

y  ( m ) 

-4854052. 280 

e Ion  ( dms  ) 

E 

276 

51 

53.99783 

1 

z  ( m ) 

4082638.525 

wlon  (dms) 

W 

83 

8 

6. 00217 

s 

ht  ( m ) 

233. 4375 

clock  off  set ( s ) 
f req  off  set ( s/s ) 

Code  ca 1 i brat i on ( m ) 
1  -  2  -.5486 

1-3  .1161 

1  -  4  -.2448 

1  -  5  -.8381 

Station  2 
Pdop  13.  1 


-.  69907204D-01 
- . 26396778D-08 

Carrier  ca 1 i brat  ion ( m ) 

-.0013 
.  0000 
-.0008 
-.0022 


| 

X 

( m ) 

593259.704 

1  at 

( dms  ) 

N 

39 

57 

48. 13617 

J 

y 

( m ) 

-4859339.368 

e  Ion 

( dms ) 

E 

276 

57 

38. 13979 

z 

( m ) 

4074952.203 

wlon  ( dms  ) 
ht  ( m ) 

W 

83  2 

129.2440 

21 . 86021 

c  lock  off  set ( s ) 
f req  off  set ( s/s ) 

Code  ca 1 i brat  ion ( m ) 
1-2  . 1555 
1  -  3  1.5778 
1  -  4  -1.0203 


.  1445  1944D-03 
. 30 199690D-09 


Carr ier 


calibration(m) 
.  0003 
-.0014 
-.0001 


» t’O 


V  ■>  V  V 


10/27/87 
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TRIMBLE  NAVIGATION,  LTD 
585  NORTH  MARY  AVENUE 
SUNNYVALE,  CALIFORNIA  94086 
U.S. A. 

PROGRAM  TRIMVEC 

GPS  RELATIVE  POSITIONING  SOLUTION 
VERSION  87. 106MB 


£>,</TTCAJ  -  IS 
■ZEWUAJ  ••*> 


File’name:  aristov.fix 
Coordinate  system  -  WGS-84 

Type  solution:  Double  difference 


Start  date/time:  1987/  3/  6 
Stop  date/time:  1987/  3/  6 

Data  available: 
station  ••  1 

sat :  9  ! . 

sat  :  1 1  ! . 

sat : 12  ! . 

sat : 13  I . 

sat :  3  ! . 


8:47:45.  day  of  year  65  tow  463665. 

10:16:60.  day  of  year  65  tow  469020. 


station : 


Ephemer  i  s 

file  used:  a:britb065. 

eph 

SATELLITE 

AODE ( hr . )  1 

HEALTH 

WEEK 

NO. 

TOW ( sec )  URA(n) 

9 

8.53 

0 

373 

463530.00 

8.0 

1 1 

11.95 

0 

373 

463500.00 

5.  7 

12 

6.26 

0 

373 

463530.00 

5.7 

13 

3.41 

0 

373 

463530.00 

5.7 

3 

4.55 

0 

373 

463560.00 

4.0 

Broadcast 

satellite  clock 

correct  ion 

values 

prn 

afO 

af  1 

af  2 

toe 

9 

. 33285003 16D-04 

.  16029844  12D- 10 

-.  2775557562D-  16 

. 4680D+06 

1 1 

. 5980860442D-04 

- . 54569682 1  ID- 1 1 

-. 2775557562D-16 

. 4680D+06 

12 

. 464 1688429D-03 

. 5570655048D- 1 1 

.  OOOOOOOOOOD+OO 

. 4680D+06 

13 

. 9786570445D-04 

. 2387423592D- 1  1 

.0000000000D+00 

. 4680D+06 

3  - 

. 3093853593D-05 

568434 1886D- 12 

.0000000000D+00 

. 4680D+06 

Message  file  for  station  1 


Station  ID:  BRITB065 
Date:  MAR. 6, 1987 


Avyvvvv 


CWi 


raw 


Pi 


4- 


Session :  065B 

Location:  BRITTON  RD.  HILLIARD 
Observer’s  name:  BRIG 
Receiver  serial  #:  4871 
Antenna  serial  #:  0117 
External  computer  serial  #:  864C 

Antenna  height:  1.3122  meters  (51.661  inches)  -  uncorrected  to  the  vertical 
BRITTON  RD.  HILLIARD 
STA.  BRIT 
SESSION  B 


Origin  of  station  1  coordinates  :  Best  C/A  code  tracking  solution 
STATION  (mark)  1 


input  data  file  1  :  a : br i tb065 . dat 
antenna  height(m)  1.276 


met 

values 

used:  pressure(mb) 

1010. 0 

temperature ( deg 

C) 

20.0 

relative  humidity(%) 

50.0 

X 

( m ) 

584399.071  lat 

( dms ) 

N 

40  3 

10.65202 

y 

( m ) 

-4854054.663  elon 

( dms ) 

E 

276  51 

54.09583 

z 

( m ) 

4082640.924  ulon 

( dms  ) 

W 

83  8 

5.90417 

ht  ( m ) 

237.0324 

Message  file  for  station  2 


Station  ID:  1883b065 
Date:  03-06-87 
Session:  b065 

Location:  mon  1883  at  corner  of  central  8<mckinely 

Observer’s  name:  dip 

Receiver  serial  #:  4872 

Antenna  serial  #:  2704a00128 

External  computer  serial  #:  5501536c 

Antenna  height:  1.5420  meters  (60.709  inches)  -  uncorrected  to  the  vertical 


STATION  (mark) 


input  data  file  1 
antenna  height(m) 


1883b065 . dat 
1.512 


met 

values  used:  pressure(mb)  1010.0 

temperature ( deg  C) 

20.0 

relative  humidity(%) 

50.0 

X 

( m ) 

593200.285  lat  (dms) 

N 

39 

57 

48. 49871 

y 

( m ) 

-4859380.869  elon  (dms) 

E 

276 

57 

35. 44279 

z 

(m) 

4074995.281  ulon  (dms) 

W 

83 

2 

24.55721 

ht  ( m ) 

182. 

9681 

Vector  1 

originates  at  station  1  ends  at 

station  2 

dz (01 ) 


i.  Vector  covariance  matrix  (cm**2)  : 

dx(01)  dy(01) 

dx(01)  .0003795 

dy(01)  -.0008262  .0074935 

dz(Ol)  .0000941  -.0021332 

Vector  correlation  matrix  : 

dx(01)  dy(01) 


0014925 
dz ( 0  1 ) 


dx(01)  1.0000000 

dy(01)  -.4898924  1.0000000 

dz(01)  .1250329  -.6378768  1.0000000 


STATION  1  TO  STATION  2 
Kdop  is  .080179 

slope  distance  (m)  12817.3983  sigma  (cm)  .037 

normal  section  azimuth  (dms)  140  47  47.40 
vertical  angle  (dms)  0-17  57.48 


east(m)  north(m)  up(m) 

Delta  lat(dms)  0  -5  22.15332 

Delta  lon(dms)  0  5  41.34697 

Delta  ht ( n )  -54.0642 

correlations  for  baseline  1: 


8101. 467 


-9932.  142 


-66. 955 


dx 

dy 

dz 

trop 

bias  1 

bias  2 

bias  3 

bias  4 

dx 

1 . 0000000 

dy 

-.  4898924 

1.0000000 

dz 

.  1250329 

-.  6378768 

1 . OOOOOOO 

trop 

.  0000000 

.  OOOOOOO 

. OOOOOOO 

1.  OOOOOOO 

bias  1 

.  OOOOOOO 

. OOOOOOO 

. OOOOOOO 

.  OOOOOOO 

1.  OOOOOOO 

bias  2 

.  0000000 

.  OOOOOOO 

. OOOOOOO 

. OOOOOOO 

.  OOOOOOO 

1.  OOOOOOO 

bias  3 

. OOOOOOO 
1.0000000 

. OOOOOOO 

. OOOOOOO 

.  OOOOOOO 

.  OOOOOOO 

.  OOOOOOO 

bias  4 

. OOOOOOO 

. OOOOOOO 

. OOOOOOO 

.  OOOOOOO 

.  OOOOOOO 

.  OOOOOOO 

. OOOOOOO 

1.0000000 

dx  ( m ) 
dy  (a) 
dz  ( a ) 
trop  (%) 
bias  1  (cycle) 
bias  2  (cycle) 
bias  3  (cycle) 
bias  4  (cycle) 


Solution 
8801. 214 
-5326.207 
-7645. 643 
.000 
.  000 
.000 
.  000 
.000 


Sigma 
.  000 
.  001 
.  000 
.  000 
.  000 
.000 
.  000 
.  000 


Results  of  integer  bias  search: 

1.76043  8.66400 

prn  11  0 


10. 90994 


Ratio  sum-o f-square s ( 2 )  to  sum-o f-squares (  1 ) 


4.92 


Interval  between  epochs  (sec)  15 
Epoch  increment  1 

Number  of  measurements  used  in  solution 
Number  of  measurements  rejected 
RMS  (cycles)  .024 

Elevation  mask  (deg)  15.0 

Edit  multiplier  3.5 

Modified  Hopfield  troposphere  model  used 

Best  tracking  C/A  code  positions 
Station  1 
Pdop  6.  1 


x  ( m ) 

584399.  188 

1  at  ( dms ) 

N  40  3 

10.65202 

y  ( m ) 

-4854055.632 

elon  (dms) 

E  276  51 

54. 09583 

z  ( m ) 

4082641.745 

wlon  (dms) 

W  83  8 

5.90417 

ht  ( m ) 

238.3084 

c lock  off  set ( s  ) 

•.  68906956D-01 

freq  of  f  set ( s/s  ) 

. 56754055D-09 

Code  cal ibrat ion ( m ) 

Carrier  ca 1 ibrat ion ( m ) 

1  -  2 

-.4253 

-.0019 

1  -  3 

.  1072 

-.0005 

1  -  4 

-.0536 

-.0012 

1  -  5 

-.6754 

-.0033 

Station 

2 

Pdop 

6.  1 

x  (m) 

593201.878 

lat  ( dms ) 

N  39  57 

48.32256 

y  ( m ) 

-4859412. 488 

elon  (dms) 

E  276  57 

35. 34799 

z  ( m ) 

4075014.655 

wlon  (dms) 

W  83  2 

24.65201 

ht  ( m ) 

219.6162 

clock  offset(s)  -.  69 159230D-0 1 

freq  offset(s/s)  - . 19 168765D-09 


Code  cal ibrat i on ( m ) 
1-2  .  1447 
1  -  3  -.2377 
1-4  -.  1537 
1  -  5  -.0929 


Carrier  ca  1  i  brat'ion  (  m  ) 

.  0003 
-. 0012 
-.0000 
-.0013 


1060 

129 


